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ABSTRACT 

Magnetohydrodynamics in divergence form describes a hyperbolic system of 
covariant and constraint-free equations. It comprises a linear combination of 
an algebraic constraint and Faraday's equations. Here, we study the problem 
of well-posedness, and identify a preferred linear combination in this divergence 
formulation. The limit of weak magnetic fields shows the slow magnetosonic and 
Alfven waves to bifurcate from the contact discontinuity (entropy waves), while 
the fast magnetosonic wave is a regular perturbation of the hydrodynamical sound 
speed. These results are further reported as a starting point for characteristic 
based shock capturing schemes for simulations with ultra-relativistic shocks in 
magnetized relativistic fluids. 

Subject headings :MHD, divergence form, characteristics 



1. Introduction 

Highly relativistic astrophysical fluids have been observed as highly energetic outflows, 
e.g.: jets in active galactic nuclei, including a few optical radio-jets such as 3C273 (Pearson 
et al. 1981, Thomson et al. 1993, Bhacall et al. 1995), 3C346 (Dey & van Breugel 1994), 
M87 (Biretta et al. 1995) and PKS 1229-21 (Le Brun et al. 1996), microquasars in our 
galaxy (Hjellming & Rupen 1995; Mirabel & Rodriguez 1995; Levinson & Blandford 1996), 
pulsar winds (Kennel & Coroniti 1984), and fireballs in recent models of 7— ray bursts (Rees 
& Meszaros 1995). These flows are generally time-dependent, or have been produced in a 
strongly time-variable episode, and hence are relativistically shocked fluid flows. In most 
cases, shocks are responsible for brightest emission features at the highest energies. 
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The evolution of strongly magnetized flows can be markedly different from 
unmagnetized flows. This is already apparent from in small amplitude wave-motion in ideal 
magnetohydrodynamics compared with hydrodynamics, and their distinct shock structures. 
The nonlinear development of large scale morphology of strongly magnetized jets can result 
in features such as the formation of a nose cone (Clarke et al. 1986), which is absent in 
hydrodynamical evolution. Of particular interest is the role of magnetic fiels in the large 
scale, three-dimensional stability of jets and their knotted structures. 

Time-dependent simulations may provide the link between the observed emission 
features and the internal structure such as magnetized field distribution, and boundary 
conditions at the source. It is hoped that simulations ultimately provide constraints on 
the flow parameters, perhaps also derived from stability criteria. Higher dimensional 
simulations of jets are performed by a number of groups in the approximation of relativistic 
hydrodynamics (van Putten 1993b; Duncan & Hughes 1994; Marti et al. 1995, Gomez et 
al. 1997) and relativistic magnetohydrodynamics (van Putten 1994ab, 1996; Nishikawa et 
al. 1997; Koide et al. 1996, 1998). 

The earliest approach for time-dependent simulations on shocked relativistic 
magnet ohydro dynamic flows with dynamically significant magnetic fields uses the equations 
of magnetohydrodynamics (MHD) in divergence form (van Putten 1991, 1993a). The 
divergence technique obtains hyperbolic systems from partial differential-algebraic systems 
of equations, and applies more generally to the case of Yang-Mills magnetohydrodynamics 
in SU(N) (van Putten 1994cd, Choquet-Bruhat 1994ab), and general relativity (van Putten 
& Eardley 1996). A linear smoothing method has been used as a shock capturing scheme 
for this formulation (van Putten 1993a, 1994ab, 1995). Both one- and two-dimensional 
simulations on astrophysical jets are performed (van Putten 1993b, 1996, Levinson & 
van Putten 1997). This method is accurate and stable, and generally performs well for 
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relativistic shocked fluid flow with up to moderately strong shock strengths (van Putten 
1993), and preserves divergence free magnetic fields to within machine round-off error (van 
Putten 1995). A smoothing method, therefore, is appropriate for simulations on the large 
scale morphology of astrophysical jets. 

Advanced shock capturing schemes are commonly based on characteristics, however, 
such as Roe's method (1981) and its extensions. These methods are generally more stable 
than smoothing methods for flows with ultra-relativistic shocks, such as in calculations 
of fire-balls for 7— ray bursts (Rees & Meszaros 1994; Wen et al. 1997). It is therefore of 
interest to explore applications of these shock capturing schemes to relativistic MHD. Here, 
we describe a first step in this direction is given by studying the computational stability of 
normalized right nullvectors (i.e., the right eigenvectors) of the characteristic matrix. 

The divergence technique incorporates a constraint c = into a divergence equation of 
a two-form, V a u a b = 0, as in Faraday's equation, through the linear combination 

V a (uj ah + Xg ah c) = 0, A^O. (1) 

In the context of an Cauchy problem, (1) conserves c = in the future domain of dependence 
of the initial hypersurface with physical Cauchy data (van Putten 1991). 

In this paper, we identify a preferred linear combination in (1), i.e.: a choice of A 
and overall sign of (1)) in its application to the equations of ideal MHD. This follows 
from two separate analysis: a derivation of the right nullvectors of the characteristic 
matrix and well-posedness. Somewhat remarkably, both analysis agree in their preferred 
linear combinations. This suggests to consider this preferred linear combination in future 
applications in characteristic based methods to MHD in divergence form. 

The problem of linearized perturbations in relativistic MHD has been considered 
previously by Anile (1989) and that of Alfven waves by Kommissarov (1997). The 



- 5 - 



well-posedness proof uses an extension to the Friedrichs-Lax symmetrization procedure 
from earlier work on Yang-Mills magnetohydrodynamics (van Putten, 1994cd). 

Section 2 describes non- uniqueness in the original version of the divergence technique. 
In Section 3, a new derivation of the right nullvectors is given. Section 4 briefly summarizes 
well-posedness obtained by embedding of physical solutions in a symmetric hyperbolic 
system of equations. 



2. MHD in divergence form 

Ideal MHD describes an inviscid, perfectly conductive plasma in a single fluid 
description with velocity four- vector, u b {u c u c = — 1). It is given by energy-momentum 
conservation, V ' a T ab = 0, where T ab is the stress-energy tensor of both the fluid and the 
electromagnetic field, Faraday's equations, V a (w' a /^) = subject to the algebraic constraint 
u c h c = 0, and conservation, V a (ru a ) = 0, of baryon number, r. For a polytropic equation of 
state with polytropic index 7, we have T ab = (r + ^pjj + h 2 )u a u b + (P + h 2 /2)g ab — h a h b , 
P is the hydrostatic pressure and g ab is the metric tensor. The theory of relativistic 
magnetohydrodynamics is contained in the conservation laws of energy-momentum, 
V a T ab = 0, and baryon number, V a (rw a ) = 0, together with Faraday's equations and a 
constraint, 

V a (h [a u b] ) = 0, u c h c = 0. (2) 

The divergence technique considers a constraint-free formulation by taking a linear 
combination 

V a {h [a u b] + \u c h c ) = 0. (3) 

Provided A 7^ 0, (3) preserves u c h c = during dynamical evolution in response to physical 
initial data (van Putten 1991), and no constraint violating wave-motion occurs. 
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Algebraically, the linear combination (3) establishes a rank-one update to its Jacobian, 
and hence of that of the full equations of MHD. Clearly, symmetry conditions of the 
Jacobian may enter a particular choice of A. Below, we consider the choice 

A = 1, (4) 

so that 

V a T ah = 0, 
-V a (h^vfi + g ab u c h c ) = 0, 

(5) 

V a {ru a ) = 0, 
V a (eV + l)) = 0, 

where £ is any time-like vector field and U = (u b , h b , r, P). The minus sign in front of the 
present linear combination is chosen also in regards to the structure of the Jacobian of (5). 
This will be made explicit below. 

Upon expansion, (5) obtains the system 

A a d a U + • • • = 0, (6) 

where the matrices A a % = A°^(U) = are 10 by 10, and the dots refer coupling terms 
to the Christoffel symbols. The infinitesimal wave-structure is given by characteristic 
wave-fronts at given U (since the A a are coordinate independent). The simple wave ansatz 
U = U(4>) obtains 

A a d a (j)U' + • • • = 0. (7) 

The wave-fronts are characteristic surfaces, whenever the matrix A a d a (fi is singular. 
The directions v a = d a (j) then are the normals to these surfaces. The small amplitude 
perturbations in these simple waves are given by the right nullvectors of A a u a . Stated 
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differently, the small amplitude perturbations are right eigenvectors, R, of (A 1 ) 1 A x u x , 
when the wave moves along the x— direction, following 

((A t )' 1 A x - v) R = 0, (8) 

where v is the velocity of propagation. 

The divergence technique provides an embedding of the theory of ideal MHD in a 
system of ten equations. Physical initial data are properly propagated by it, without 
exiting non-physical wave-modes. The physical waves (entropy waves, Alfven and 
magnetohydrodynamic waves) are all contained within the light cone. Here, adding 

g ab u c 

h c (or a multiple thereof) to Faraday's equations provides a rank-one update to the 
characteristic matrix A c v c . On the light cone, however, v 2 = 0, and this linear combination 
no longer regularizes of the characteristic determinant. (This results from insisting on 
covariance in the divergence formulation.) Kommissarov (1997) attempts to discusses MHD 
in divergence form outside the context of the initial value problem with physical initial 
data, and hence erroneously concludes the presence of non-physical wave-modes. 



3. The characteristic matrix 



We have 



dF 



aA 



Q F aA 



QJJB - QJJB 



with p = r + — + h 2 = rf + h 2 , we have 



F cA u c = < 



p{u c v c )u a + (P + h 2 /2)u a - (h c u c )h a , 
-{{h c v L )u a - (u c v c )h a + u a u c h c }, 

(^ c )(u 2 + l). 



(9) 



(10) 
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The system of 10x10 equations for U B = (u b , h b , r, P) can be reduced to 8x8 in the variables 
V B = (v s ,h b ,r) by expressing u h in terms of the spatial three- velocity u b = T(l,v s ), 
T = 1/Vl — v 2 (l, v s ), s = 1,2,3, noting that linearized wave-motion conserves entropy, so 
that dP = 7^dr. In V B , the equation of energy conservation, V a T at = and the last 
equation of (5) are automatically satisfied, whence they can be ignored. In what follows, A a 
shall denote the resulting 8x8 matrix, obtained from the original 10x10 matrix by deletion 
of the first and last row, addition of the last column (multiplied by -y-) to the one-but last 
column (associated with r), followed by deletion of the first and last columns. 

The linearized wave-structure is given by the characteristic problem 

A c u c z = (11) 

for the right null-vectors z = U' . Without loss of generality, (11) can be studied in a 
co-moving frame, in which u b = (1, 0, 0, 0). In this event, T = 1 and ^7 = 0. Furthermore, 
the x-axis of the local coordinate system can be aligned with the magnetic field, so that 
h b = (0, if, 0,0). Given the two orientations u s and h b , the wave-structure is rotationally 
symmetric bout the x-axis, and hence u y and u z act symmetrically as ^Jv 2 + v 2 z \ we will put 
v z = 0. For A c v c , we have 



pv x -vxtL -Hv 2 -Hv 3 ' r 2 

p V1 Hv 3 -Hv 2 ~' P r " 3 

pvi -Hv 2 

v t H — v\ -v-2 -v 3 

(12) 

— Hv 2 Hv 3 v 2 vi 

-Hv 3 -Hv 2 v 3 i>i 

-Hv 2 v t 

rv 2 rv 3 v\ 



Note that the lower diagonal block is v\ times the 4x4 indentity matrix. This results 
from the sign choice in the given combination of Faraday's equations and the constraint 
in (5) and (10). Furthermore, notice that the third and seventh rows and columns act 
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independently to give rise to the Alfven waves. The remaining waves are described by the 
reduced problem 



(A c u c )'z' = 0, 



(13) 



where (A c z/ C )' is obtained from A c u c by deleting the third and seventh rows and columns, 
thereby obtaining a problem in the 6-dimensional variable z' . Introducing 



y 

[11) takes the form of a coupled system of 3x3 equations 

V\Zx + Xy = 0, Yx + Viy = 0, 



in which 



p -H 
p 
H -1 





-Hv 2 


-Hu 3 


7 Pv'i 
r 


X = 


Hu 3 


-Hu 2 


jPl>3 
r 




-v 2 


-u 3 







-Hu 2 


Hv 3 


v 2 


Y = 


-Hv 3 


-Ev 2 


vz 




rv 2 


rv 3 






This obtains a single 3x3 eigenvalue problem in x, given by 

XYx = v\Zx & Z~ 1 XYx = ufx. 



(14) 



(15) 



(16) 



(17) 
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Here, Z 1 XY — v\ is given by the matrix 



Wi,i 


w 1:2 


W 2 ,i 




rf 


H-, Pv 2 v 3 
rf 






2,2 2 
v 2 + v 3 - v\ 



where the upper diagonal 2x2 matrix W is given by 

W = Wij. (19) 

The two zeros in the third column of (18) result from A = 1. Upon substitution 
z/| = v 2 + v\ — u%, the determinant assumes the covariant expression 

pdetW= (r/- 7 P)(w c z/ c ) 4 



-(h 2 + 1 P)u 2 (u c u c ) 2 + ^(h c u c ) 2 u 2 . 
Alfven waves. The eigenvalues for the Alfven waves are given by 

\h c u r \ 



(20) 



(21) 



with null-vector 

z = (0,0, #^,0,0,0,^, Of, 
associated with Alfven waves; covariantly, 

U A = (v a ,±^v a ,0,0) T , 

where v a may be taken to be 

H(Q, 0, i/ 4 , -i/ 3 ) = e abcd u h h c u d = v a . 
Thus, the Alfven wave is transversal in which h 2 is conserved (#/i fe is orthogonal to h b ). 



(22) 



(23) 



(24) 



Magnetohydrodynamic waves. The eigenvalues for the magnetohydrodynamic waves 
are given by the roots of the characteristic determinant (20). Writing 



n h = v h + (u c v c )u\ 



(25) 



we have v 2 = —t 2 + n 2 , t — u c u c , n 2 = n c n c . Let a — and (3 — ^p. Then 

(h c u c ) 2 (3 (h c n c ) 2 _P cqz 2j 
rfn 2 a h 2 n 2 a 

Consequently, (20) becomes 

(a-iy- (l + P)v 2 (l-v 2 ) 

+(3a~ 1 cos 2 0(1 - v 2 ) = 0, 

where v 2 = (27) has real solutions v for any given n b , whenever 

(a + fi)v A -{1+13 + fia-^v 2 + (3a' 1 = 
has real solutions v. But (28) has discriminant 

D = (a + (3-a(3) 2 > 0. 
Weak magnetic fields are described by small (3 expansions as follows. 

Proposition 3.1. Fast magnetosonic waves are a regular perturbation of sound waves in 
pure hydrodynamics, while the Alfven and slow magnetosonic waves bifurcate from entropy 
waves (contact discontinuities), whose propagation velocities satisfy 

vj/v 2 -l + ^sin^ + OOn 

v\/v 2 ~ /3cos 2 0[l-/ta -1 + 0(/3 2 )], (30) 

v 2 /v 2 ~ (3 cos 2 0[1 - (3(1 - 2=1 cos 2 0) + 0(/3 2 )], 
where v\ = a^ 1 is the square of the hydrodynamical velocity, and which obey the inequalities 

v 2 s <v\< v 2 . (31) 

Inequalities (31) remain valid for general (3 (e.g. Bazer & Ericson 1959; Lichnerowicz 1967; 
Anile 1989). 



(26) 
(27) 

(28) 
(29) 
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4. Right nullvectors 



Inspection of (19), together with (15), shows the null- vector 



/ 



\ 



(32) 



^1^2 ^3 



HV\V 2 v\ 

Ev 2 {v 2 -av\) 
-Ev 2 vj,{v\ - oivl) 


Of course, (32) can be stated covariantly by noting that H 2 = h 2 , Hu 2 = h c u c , v x = u c u c , 

H 2 {v\ - av\) = (h c u c ) 2 - ah 2 (u c u c ) 2 = h 2 k u (33) 

and introducing 

if (0, v\ + is 2 , -is 2 is 3 , -v 2 Va) T = e ahcd u b v c v d = w a . 
Since —arv\v\ is a scalar, u 3 is to be treated as 

H 2 {v\ + v\) = h 2 n 2 - (h c u c ) 2 = h 2 k 2 , 
were n a = u a + {u c u c )u a . Note that 

ki = n 2 (cos 2 — at; 2 ), k 2 = n 2 sin 2 <j>, 

where v = v s ,Vf. Clearly, z is formed from 

Su b = -t{k in b - (k 2 + k 1 )(h c n c )h b ) 

Sh b = hw b + k 2 t{h c n c )u b , 

5r = —ark 2 t 2 , 

SP =-rfk 2 t 2 , 



(34) 



(35) 



(36) 



(37) 



where h = h b /\h\, and 

v a = e abcd u h h c v d , w a = e abcd u b v c v d . 
We thus have the following. 



(38) 



Proposition 3.2. Given a unit vector n b orthogonal to u b , and a root v b = n b +vu b , v = u c v c 
of (28), the right nullvectors for the hydrodynamical waves of (11), U A = (5u b ,5h b ,5r,SP), 
are 



Su b = 



sin 2 n b — (1 — av 2 ){n b — cos0 h b ) 



Sh b = \h\ [(cos 2 — av 2 )w b + vsin 2 0cos0 u b ] , 



(39) 



5r = —v 2 ar sin 2 0, 



5P = —v 2 rf sin 2 (j>. 

where w b = w b /\h\ . 



Anile (1989) gives a different form of these right nullvectors. By Proposition 3.1, our weak 
magnetic field limits show that 

cos 2 - avj < (40) 

for fast magnetosonic waves, while 

cos 2 - av 2 > (41) 

for slow magnetosonic waves. Inspection of (34) shows that therefore the tangential 
component of the magnetic field is strengthened in fast magnetosonic waves, while it 
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is weakened in slow magnetosonic waves. This distinguishing aspect of fast and slow 
magnetosonic waves was first noted by Bazer & Ericson (1959) in their analysis of shocks 
in nonrelativistic MHD. 

The limit of small /3 is of particular interest to computation. For example, in 
various settings a magnetized fluid streams into a nearly unmagnetized environment. A 
characteristics based scheme must therefore reliable treat a large dynamic range in (3. 
Clearly, a full set of nullvectors (including those of contact discontinuities) obtains for 
nonzero (3. However, the behavior of these nullvectors is somewhat nontrivial as (3 becomes 
small. In what follows, we consider the small f3 limit, in the sense of small \h\/y/jP, while 
keeping the direction h constant. In this limit, 

\-av 2 ~ -/3^sin 2 + O(/3 2 ), 

y 1 (42) 

1 - av 2 ~ 1 + 0((3) 
for the fast and slow magnetosonic speeds, respectively. 



Corollary 4.1. In the limit of low magnetic field strength, the fast magnetosonic waves 
are described by the right nullvectors 

Su b =v f n b + (3^(n b -cos <f>h b )v f + 0((3 2 ), 
Sh b = \h\(—w b + Vf cos0 u b ) + (3^-w b + 0(/3 2 ), 
Sr = —v 2 ar, 
SP =-v}rf, 
and the slow magnetosonic waves by 

5u b = cos0(/i & - cos0n 6 ) + 0{(3), 

Sh b = v^P(cos0 w b + v s sm 2 (f) u b ) + 0((3), 

(44) 

Sr = —v s ar sin 0, 
SP = —v s arf sin 2 0. 
The small (3 limit of the nullvectors can now be normalized. 
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4.1. Bifurcations from entropy waves 

The behavior of the nullvectors in the limit of weak magnetic fields can be derived 
from (23) and Corollary 4.1. To this end, note that 

v a = \h\v a = sin0 \h\v a , (45) 

where v c v c = 1, and <fi denotes the angle between rf and h c , 

n b = cos0 h b + sm<f) y\ (46) 

y c u c = h c y c = 0,y c y c = 1 (n b is normalized to be unit, as in the assumptions of Proposition 
3.2). It follows that the Alfven nullvectors may be normalized to 

5U A = (v a ,± v ^v a , 0,0). (47) 

In the limit of vanishingly small f3, the pair of slow magnetosonic waves collapse to the 
single normalized nullvector 



5U A =(y b , vW,0,0). (48) 

Note that y c v c = 0, so that (47) and (48) are independent. Division by sin0 thus provides 
a normalization of the original expressions (23) and (44). 

The nullvector associated with entropy waves (u c u c = 0) is 

SU A = (0, 0, 5r, 0) (49) 

if h c u c 7^ 0, and 

(0,5h c ,5r,5P), (5u c , 0,0,0), (50) 

if h c v c = 0, subject to 

5P + h c 5h c = 0, u c 5h c = 0, v c bu c = 0. (51) 
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The second case refers to transverse MHD for which continuity must hold of total pressure, 
zero orthogonal magnetic field and transverse velocity. Note that transverse MHD has two 
nullvectors, and corresponds to the case of pure hydrodynamics. With the exception of 
transverse MHD, therefore, the contact discontinuity provides one nullvector. 

Transverse MHD or pure hydrodynamics allows for shear along contact discontinuities, 
which is responsible for the two independent nullvectors. Whenever magnetic field lines 
cross a contact discontinuity, however, the persistent coupling to the magnetic field lines in 
ideal MHD prohibits shear. In ideal MHD, the response to the original two-dimensional 
degree of freedom in shear is two new wave-modes. These wave modes are the Alfven 
wave and the slow magnetosonic wave. These two wave- modes are indeed different, as (47) 
and (48) show. The Alfven and slow magnetosonic wave may be regarded as one pair, 
bifurcating from the contact discontinuity. This has been illustrated in Fig. 6 of van Putten 
(1993a). Indeed, the limit of vanishing f3 recovers the two shear modes from the independent 
Alfven and slow magnetosonic waves. Of course, the Alfven wave is purely rotational, 
while the slow magnetosonic wave is slightly helical, including a longitudinal variation of 
±v s sin 2 = ±/3 sin 2 0cos0. The fast magnetosonic wave remains a regular perturbation of 
the ordinary sound wave. 

The weak magnetic field limit thus obtains two nullvectors from the fast magnetosonic 
waves, two from the Alfven waves, one from the slow magnetosonic waves and generally one 
from the contact discontinuity, a total of six. This leaves an apparent degeneracy of one. 

The degeneracy stems from the neighboring to order v s of the two nullvectors of the 
slow magnetosonic waves. This would suggest ill-posedness to this order in projections. 
However, characteristic based methods consider the product of the projections on the 
nullvectors and the associated eigenvectors. In the present case, therefore, the order of 
the degeneracy is precisely cancelled by multiplication with the eigenvalue v s , which is 
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computationally stable. The limit of arbitrarily small (3 in the application of characteristic 
based methods is computationally well-posed. 

5. Well-posedness 

The theory of ideal relativistic MHD was first shown to be well-posed by Friedrichs 
(1974), using the Friedrichs-Lax symmetrization procedure (1971). The problem of 
constraints was circumvented by a reduction of variables. The symmetrization procedure of 
Friedrichs and Friedrichs and Lax (1971) applies to hyperbolic systems of equations of the 
form 

V a F aB = f (52) 

which satisfy a certain convexity condition. The presence of conserved constraints, however, 
can be treated also by an extension of the Friedrichs-Lax symmetrization procedure, with 
no need for an additional reduction of variables, developed in earlier work on Yang-Mills 
magnetohydrodynamics in SU(N) (van Putten 1994cd). Once in symmetric hyperbolic form, 
well-posedness results from standard energy arguments (e.g. Fisher & Marsden 1972). The 
main arguments of symmetrization in the presence of constraints are briefly recalled here, to 
highlight the same linear combination of (5), now from the point of view of well-posedness. 



5.1. Symmetrization with constraints 

Variations 5V A of (u b , h b , r, P) can be unconstraint (with respect to all ten degrees 
of freedom), and constraint, i.e., those obeying the constraints. For example, 5c ^ 
results from a total variation, while 5c = is a constraint variation. Symmetrization in 
the presence of constraints follows if there exists a vector field Wa which produces a total 
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derivative in the modified main dependency relation 

YI : W A 5F aA = 5z a , (53) 

and which obtains constraint positive definiteness in 

YII : 5W A 5F a % > (54) 

for some time-like vector £ a . Of course, the source terms f B must satisfy the consistency 
condition 

W A f A = (55) 

whenever the constraints are satisfied. Allowing a possible nonzero total derivative in YI 
defines an extension (van Putten 1994cd) to the Friedrichs-Lax (1971) symmetrization 
procedure. 

Differentiation by V c of the unconstraint identity YI obtains 

V D + WAd2Fa \ V o = d2z y V d m 
dV c dV D dV c dV D " dV c dV D ' y ' 

This establishes symmetry of the matrices 

_dW A dF« A 

Also, 

6V c A a CD tJV D =r§)(i?^)=Wt>0 (58) 

for all constraint variations 5V A . Of course, given V A , the constraint variations 5V A define 
a linear subspace V of dimension N — m, where m is the number of constraints c = 0, each 
giving rise to 

= So = JtL5V\ (59) 
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We have the following construction (van Putten, 1994cd). 

Lemma 5.1 Given a real- symmetric AeC(R n , R n ) which is positive definite on a linear 
subspace V C R™ ; there exists a real-symmetric, positive definite A*e£(R n , R n ) such that 



This may be seen as follows. Consider A* = A + fix T x, where a; is a unit element 
from V L . Then A* is symmetric positive definite on V — {z — y + \x\yeV, XeR} : 
z T A T z > c'\\z\\ 2 = c'(\\y\\ 2 + X 2 \\x\\ 2 ) with d > upon choosing fi> M, where M — \ \A\\ 
denotes the norm of A. This construction may be repeated until V 1 - is exhausted, leaving 
A* symmetric positive definite on R n as an embedding of A on V. 

The real-symmetric matrix A^j D ^ a is positive definite on the subspace of constraint 
variations V; let {A^ D ^ a )* be the positive definite, symmetric matrix obtained from the 
Lemma. It follows that solutions to (52) (and its constraints) satisfy the symmetric positive 
definite system of equations 



A*y = Ay (yeV). 



(60) 



(A aAB yaev c )v A + 



A aAB {V,:) a V A = f B , 



(61) 



where 



v a = -ar v c ) + (v s ) a . 

It remains to show that ideal MHD satisfies properties YI and YII. 



(62) 
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5.2. Symmetrization of hydrodynamics 

Relativistic hydrodynamics has been shown to be symmetrizable by Friedrichs (1974), 
Ruggeri & Strumia (1981), and Anile (1989). This uses the equations in the form 



V a (rfu a u b + Pg ah ) =0, 

V a (ru a ) =0, (63) 
X7 a (rSu a ) =0 

away from entropy generating shocks. Then = (u a , f — TS, T) and Vq = (v a , T, /) with 
a reduction of variables on the velocity four- vector by u b = T(l,v a ), where T is the Lorentz 
factor. With Ff A denoting the fluid dynamical equations V a Tf = 0, Tf = rfu a u b + Pg ab 
with / the specific enthalpy, and V a (ru a ) = 0, it has been shown that (Ruggeri & Strumia 
1981; Anile 1989) 

W{5Ff = 0, Q f = 5W A 5F?%>0 (64) 

provided that the free enthalpy G(T, P) = f — TS — 1 is concave, and the sound velocity 
is less than the speed of light. Under this conditions, the hydrodynamical equations alone, 
therefore, satisfy YI and YII, and in fact the original Friedrichs-Lax conditions CI and CII 
of Friedrichs & Lax (1971), so that they satisfy a symmetric hyperbolic system of equations. 

5.3. Symmetrization of ideal MHD 

In what follows, we set 

Uab = h a u b - u a h b + g ab u c h c , 

(65) 

= h 2 u a u b + \h 2 g ab - h a h b . 
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(66) 



We then have the expansions 

utSTg = u b (h 2 u a 5u b + h 2 u b 5u a + 2u a u b h c 5h c 
+g ab h c 5h c - h a 5h b - h b 5h a ) 
= -h 2 5u a - u a (h c 5h c ) - h a (u c 5h c ) 
-c5h a , 

h b 5uj ab = h b (h a 5u b + u b 5h a - h b 5u a -u a 5h b 
+g ab 5c) 

= h a (h c 5u c ) + c5h a - h 2 5u a - u a (h c 5h c ) 
+h a 5c. 

We hereby arrive at the identity 

u h 8T% - h b 5uj ab = 5z a , (67) 

where z a = —2h a c. The total derivative in (67) follows by the unique linear combination 
h a u b - h b u a + g ab c, as in (5). With W A = {u a , h a , f - TS, S) and F aA given by (5) 



CO 



ab 



[rewritten according to (63)], it follows that 

W A 5(Ff + F^ A ) = 5z a . (68) 

A similar calculation (van Putten 1994cd) shows that quadratic of constraint variations Q m 
given by 

5u b 5T^ a - 5h b 5cu a % = (u^ c )[h 2 (5u) 2 + (5h) 2 ] 

(69) 

+ 2[{i c 5u c ){h c 5h c ) - m c )(5u c 5h c )} 
is positive definite (for 5h a ^ 0). Therefore, the sum 

Q = 5W A 5F aA £ a = Q f + Q m (70) 

is constraint positive definite, whenever Qf is such (with respect to the fluid dynamical 
variables). It follows that both YI and YII are satisfied (with W A = (u a , h a , f — TS, S) and 
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Va = (v a , h a , T, /)), and hence physical solutions to (5) satisfy the symmetric hyperbolic 
system (61) with f B = 0. 
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